A DEPARTURE FROM PREDICTION: ELECTROWEAK PHYSICS AT NUTEV 
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Abstract 

The NuTeV experiment has performed precision measurements of the ratio of neutral- 
current to charged-current cross-sections in high rate, high energy neutrino and anti- 
neutrino beams on a dense, primarily steel, target. The separate neutrino and anti-neutrino 
beams, high statistics, and improved control of other experimental systematics, allow the 
determination of electroweak parameters with significantly greater precision than past vN 
scattering experiments. Our null hypothesis test of the standard model prediction mea- 
sures sin 2 #^ n_shell) = 0.2277±0.0013(stat) ±0.0009(syst), a value which is 3.0a above 
the prediction. We discuss possible explanations for and implications of this discrepancy. 



1 Introduction and Motivation 



Neutrino scattering played a key role in establishing the structure of the Standard Model 
of electroweak unification, and it continues to be one of the most precise probes of the 
weak neutral current available experimentally today. With the availability of copious data 
from the production and decay of on-shell Z and W bosons for comparison, contempo- 
rary neutrino scattering measurements serve to validate the theory over many orders of 
magnitude in momentum transfer and provide one of the most precise tests of the weak 
couplings of neutrinos. In addition, precise measurements of weak interactions far from 
the boson poles are inherently sensitive to processes beyond our current knowledge, in- 
cluding possible contributions from leptoquark and Z' exchange ^) and new properties of 
neutrinos themselves ^) . 

The Lagrangian for weak neutral current u-q scattering can be written as 

c = -5^.(1 -yy) 

x (497„(1 - 7 5 )« + 4977.(1 + 7 5 )«) , (D 

where deviations from p — 1 describe non-standard sources of SU(2) breaking, and e\ R 
are the chiral quark couplings 1 For the weak charged current, e| = J^ ak and e q R = 0, but 
for the neutral current e q L and e R each contain an additional term, — Qsin 2 9w, where Q is 
the quark's electric charge in units of e. 

The ratio of neutral current to charged current cross-sections for either v or V 
scattering from isoscalar targets of u and d quarks can be written as ^) 

gW5 «Vn~Vx) =(9! + r (- m (2) 

a(V N -> £-(+)X) 

where 

= a(vN -> e+X) 1 
r ~ a{uN^£-X) ~ 2' PJ 

and (7^ ^ = (e'l R ) 2 + {e d L R } 2 . Many corrections to Equation 2 are required in a real 
target ^) , but those most uncertain result from the suppression of the production of charm 
quarks in the target, which is the CKM-favored final state for charged-current scattering 
from the strange sea. This uncertainty has limited the precision of previous measurements 
of electroweak parameters in neutrino-nucleon scattering ^' ^> 7) Q ne wa y to re( j uce the 

1 Note that although we use a process-independent notation here for a tree-level p, radiative corrections 
to p depend slightly on the particles involved in the weak neutral interaction. In this case, p = y/ p(")p(9). 



uncertainty on electroweak parameters is to measure the observable 

R v - ri? 57 



ir 



1 -r 



(4) 



first suggested by Paschos and Wolfenstein °> and valid under the assumption of equal 
momentum carried by the u and d valence quarks in the target. Since a uq = a uq and a vq = 
a uq , the effect of scattering from sea quarks, which are symmetric under the exchange 
q <-> q, cancels in the difference of neutrino and anti-neutrino cross-sections. Therefore, 
the suppressed scattering from the strange sea does not cause large uncertainties in R~ . 
R~ is more difficult to measure than R v , primarily because the neutral current scatterings 
of v and V yield identical observed final states which can only be distinguished through a 
priori knowledge of the initial state neutrino. 

The experimental details and theoretical treatment of cross-sections in the NuTeV 
electroweak measurement are described in detail elsewhere ^) . in brief, we measure the 
experimental ratio of neutral current to charged current candidates in both a neutrino and 
anti-neutrino beam. A Monte Carlo simulation is used to express these experimental ra- 
tios in terms of fundamental electroweak parameters. This procedure implicitly corrects 
for details of the neutrino cross-sections and experimental backgrounds. For the measure- 
ment of sin 2 9 W , the sensitivity arises in the v beam, and the measurement in the V beam is 
the control sample for systematic uncertainties, as suggested in the Paschos-Wolfenstein 
R~ of Eqn. 4. For simulataneous fits to two electroweak parameters, e.g., sin 2 W and p or 
left and right handed couplings, this redundant control of systematics cannot be realized. 

2 Result 

As a test of the electroweak predictions for neutrino nucleon scattering, NuTeV performs 
a single-parameter fit to sin 2 9 W with all other parameters assumed to have their standard 
values, e.g., standard electroweak radiative corrections with p = 1- This fit determines 

sin 2 0^ n-shell) = 0.22773 ± 0.00135(stat.) ± 0.00093(syst.) 

M 2 — (175 GeV) 2 
- ° °0022 x ( , 

+ 0.00032 x Ml gS£). (5) 

The small dependences in M top and M H igg S result from radiative corrections as deter- 
mined from code supplied by Bardin / and from V6.34 of ZFITTER 

10) 

; however, it 



should be noted that these effects are small given existing constraints on the top and 
Higgs masses A fit to the precision electroweak data, excluding neutrino measure- 
ments, predicts a value of 0.2227 ± 0.00037 1 1 > 12 ) , approximately 3a from the NuTeV 
measurement. In the on-shell scheme, sin 2 9 W = 1 — M^/M|, where M w and M z are 
the physical gauge boson masses; therefore, this result implies M w = 80.14±0.08 GeV 2 . 
The world-average of the direct measurements of M w is 80.45 ± 0.04 GeV The fact 
that the NuTeV sin 2 #^ n ~ shc11 ^ deviates so substantially from M w makes it difficult to ex- 
plain the difference between NuTeV and the standard model prediction in terms of oblique 
radiative corrections. 

Although NuTeV was primarily designed to measure sin 2 #^ n_shc11 ^ using the 
Paschos-Wolfenstein relationship, it is also possible to fit for the single parameter, p . 
As noted above, the mechanism by which the Paschos-Wolfenstein relationship reduces 
systematic uncertainties in the sin 2 9w fit is evident in the fact that R" only is sensitive 
to sin 2 9 W and thus R" essentially measures systematics common to the v and V beams. 
Because R v and R? are both sensitive to po, there is less control of theoretical systematics 
than can be achieved with the sin 2 Q w measurement, and uncertainties on p are therefore 
larger. This fit obtains 



po = 0.9942 ± 0.0013(stat.) ±0.0016(syst.) 

,M 2 o p -(175G eV) 2 
(50 GeV) 5 



+ 0.00006 x( 1 " p , en * 2 ) 



" ° 00016 X '"(fe)' <« 

Note that these two results are highly correlated; a simultaneous fit to sin 2 9 W and p 
finds: 

po = 0.99789 ± 0.00405, sin 2 9 W = 0.22647 ± 0.00311, (7) 

with a correlation coefficient of 0.850 between the two parameters. This suggests one but 
not both of sin 2 6^ n ~ shcn ^ or p may be consistent with expectations, and that NuTeV is 
unable to distinguish between these two possibilities with significant confidence. 

Finally, we have also performed a two-parameter fit in terms of the isoscalar 
combinations of effective neutral-current quark couplings (g^R) 2 = (u°^ R ) 2 + (d e ^ R ) 2 . 
The effective couplings, which describe observed experimental rates when the processes 
described by Eqn. 1 are calculated without electroweak radiative corrections, are mea- 

2 As noted above, this extraction of sin 2 6 w n shc11 ' is done assuming radiative corrections from the 
standard model as parameterized from oiem, Gf, Mz, m top and mmggs from fits to the electroweak data. 
An alternative approach, would be to fit for Mw by determining regions of m top and mjjiggs favored by 
the NuTeV data and then using those to extract the standard model prediction. 



sured at (q 2 ) ~ -20 GeV 2 . We find 3 : 



( g f f = 0.30005 ± 0.00137, {gff = 0.03076 ± 0.00110, (8) 

with a correlation coefficient of —0.017. The predicted values from Standard Model 
parameters corresponding to the electroweak fit described earlier ^> ^) are {gf 1 ) 2 = 
0.3042 and (<?jf ) 2 = 0.0301. Note that due to the asymmetry between the strange and 
charm seas and to the slight excess of neutrons in our target, the NuTeV result is weakly 
affected (~ 1/30 the sensitivity of (g e ^ R ) 2 ) by the isovector combinations of couplings, 
(3 c l,r) 2 = ( u °l,r) 2 ~ (^°l,r) 2 - The results above assume standard model values for (5^ R ) 2 . 

3 Interpretation 

The NuTeV sin 2 9 W result is approximately three standard deviations from the prediction 
of the standard electroweak theory. This, by itself is surprising; however, it is not im- 
mediately apparent what the cause of this discrepancy might be. We discuss, in turn, the 
possibility that the NuTeV result is a statistical fluctuation among many precision results, 
the possibility that unexpected parton asymmetries in the NuTeV affect the result, and 
finally possibilities for non-standard physics which could be appearing in the anomalous 
NuTeV value. 

3.1 Impact on Standard Model Fits 

For fits assuming the validity of the Standard Model, it is appropriate to consider the a 
priori null hypothesis test chosen in the proposal of the NuTeV experiment, namely the 
measurement of sin 2 6>(^ n ~ shcll \ A fit to precision data 4 , including NuTeV, has been per- 
formed by the LEPEWWG ^ ) , and the contribution of each measurement to the x 2 an d 
final Higgs mass from this fit is shown in Figure 1. The global % 2 , which has significant 
contributions from NuTeV's sin 2 9 W measurement and A°p b B from LEP I, is a rather un- 
healthy 28.8 for 15 degrees of freedom. The probability of the fit x 2 being above 28.8 
is 1.7%. Without NuTeV, this probability of the resulting x 2 is a plausible 14%. This 
suggests that in the context of all the precision data, as compiled by the LEPEWWG, the 

3 Note that these coupling results are slightly different (<< la) than the value given in the published 
NuTeV result ^) due to a numerical error. 

4 It is a tautology, but still worth noting explicitly, that certain choices about what constitutes "the preci- 
sion data" must be made in order to make such a global analysis. Some choices that are made in compiling 
this data, for example, not listing the W mass for each experiment separately, decrease the number of de- 
grees of freedom without significantly decreasing the global x 2 - Other choices taken, for example choosing 
particular re-evaluations of the central value and uncertainties iJ ' i;5 ' 

16) 

, rather than others 

17) 

in 

the atomic parity violation measurement ^) G f Qfy, decrease the global x 2 . 
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Figure 1: The precision data, as compiled and fit by the LEPEWWG The global fit 
X 2 is 28.8/15. Most of the data is consistent with a low Higgs mass, except A°p B and the 
NuTeV sin 2 6 W . 

NuTeV result is still a statistical anomaly sufficient to spoil the fit if the standard model 
is assumed. 

This large x 2 is dominated by two moderately discrepant measurements, namely 
A% and the NuTeV sin 2 6 W , and if one or both are discarded arbitrarily, then the data 
is reasonably consistent with the standard model. However, the procedure of merely dis- 
carding one or both of these measurements to make the fit "work" is clearly not rigorous. 
Furthermore, the potential danger of such a procedure has been noted previously in the lit- 
erature. For example, if A°p b B were disregarded, then the most favored value of the Higgs 
mass from the fit would be well below the direct search limits 

19) 

Motivated by this large global x 2 of the precision electroweak data, we attempt 
to find other explanations for the discrepancy in the NuTeV sin 2 6\y measurement. 



3 . 2 Unexpected QCD Effects 



As noted above, corrections to Eqns. 2 and 4 are required to extract electroweak parame- 
ters from neutrino scattering on the NuTeV target. In particular, these equations assume 
targets symmetric under the exchange of u and d quarks, and that quark seas consist of 
quarks and anti-quarks with identical momentum distributions. 

The NuTeV analysis corrects for the significant asymmetry of d and u quarks 
that arises because the NuTeV target, which is primarily composed of iron, has a 5.74 ± 
0.02% fractional excess of neutrons over protons. However, this assumption is made 

under the assumption of isospin symmetry, i.e., ^u p (x) = d n ( x ), d P {x) = ^u n (x). 
This assumption, if significantly incorrect, could produce a sizable effect in the NuTeV 
extraction of sin 2 6 W ^ 21, 22, 23) 

Dropping the assumptions of symmetric heavy quark seas, isospin symmetry 
and a target symmetric in neutrons and protons, but assuming small deviations in all 
cases, the effect of these deviations on R~ is ^4) : 

5RT » -5N[^f^^f^](3A 2 u + A 2 d ) 




^(3A 2 + A 2 ) 

u p - u p -r u p — Up) 

+ fj _nln-n (2A ^ " 3(A ^ + A ' )ec) ' (9) 

Up Up-]- Up Up 

where A 2 d = {e u £ d ) 2 — (e^' d ) 2 , Qn is the total momentum carried by quark type Q in 
nucleon N, and the neutron excess, 5N = (A — 2Z)/A. e c denotes the ratio of the 
scattering cross section from the strange sea including kinematic suppression of heavy 
charm production to that without kinematic suppression. The first term is the effect of the 
neutron excess, which is accounted for in the NuTeV analysis; the second is the effect of 
isospin violation and the third is the effect of an asymmetric strange sea. 

NuTeV does not exactly measure Rr , in part because it is not possible ex- 
perimentally to measure neutral current reactions down to zero recoil energy. To pa- 
rameterize the exact effect of the symmetry violations above, we define the functional 
F[sin 2 9 W) 5; x] such that 

Asin 2 ^= f 1 F[sm 2 9 w ,5;x}5( y x) dx, (10) 
J o 

for any symmetry violation 5(x) in PDFs. All of the details of the NuTeV analysis are 
included in the numerical evaluation of the functionals shown in Figure 2. For this analy- 
sis, it can be seen that the level of isospin violation required to shift the sin 2 9 W measured 




Figure 2: The functionals describing the shift in the NuTeV sin 2 6w caused by not cor- 
recting the NuTeV analysis for isospin violating u and d valence and sea distributions or 
for (s(x)) 7^ (s(x)). The shift in sin 2 6 W is determined by convolving the asymmetric 
momentum distribution with the plotted functional. 



by NuTeV to its standard model expectation would be, e.g., D p — U n ~ 0.01 (about 
5% of D p + U n ), and that the level of asymmetry in the strange sea required would be 
S - S ~ +0.007 (about 30% of S + S). 

3.2.1 Isospin Violations 

Several recent classes of non-perturbative models predict isospin violation in the nu- 
cleon 21, 22) ear ij est estimation in the literature, a bag model calculation 20), 
predicts large valence asymmetries of opposite sign in u p — d n and d p — u n at all x, which 
would produce a shift in the NuTeV sin 2 9 W of —0.0020. However, this estimate neglects 
a number of effects, and a complete bag model calculation by Thomas et al. 21) con- 
cludes that asymmetries at very high x are larger, but the asymmetries at moderate x are 
smaller and of opposite sign at low x, thereby reducing the shift in sin 2 Q w to a negligible 
—0.0001 . Finally, the effect is also evaluated in the Meson Cloud model 22) ; and there the 
asymmetries are much smaller at all x, resulting in a modest shift in the NuTeV sin 2 6 W 
of +0.0002. 

Models aside, the NuTeV data itself cannot provide a significant independent 
constraint on this form of isospin violation. However, because PDFs extracted from neu- 
trino data (on heavy targets) are used to separate sea and valence quark distributions which 
affect observables at hadron colliders ^\ global analyses of PDFs including the possi- 
bility of isospin violation may be able to constrain this possibility experimentally. 

3.2.2 Strange Sea Asymmetry 

If the strange sea is generated by purely perturbative QCD processes, then neglecting 
electromagnetic effects, one expects (s(x)) = (s(x)). However, it has been noted that 
non-perturbative QCD effects can generate a significant momentum asymmetry between 
the strange and anti-strange seas 27, 28, 29) 

By measuring the processes v^,V N — > fi + fi~X the CCFR and NuTeV exper- 
iments constrain the difference between the momentum distributions of the strange and 
anti-strange seas. Within the NuTeV cross-section model model, this data implies a neg- 
ative asymmetry 24) ; 

S-S = -0.0027 + 0.0013, (11) 

or an asymmetry of 11 ± 6% of (S + S). Therefore, dropping the assumption of strange- 
antistrange symmetry results in an increase in the NuTeV value of sin 2 9 W , 



Asin 2 9 W = +0.0020 ± 0.0009. 



(12) 



The initial NuTeV measurement, which assumes (s(x)) = (s(x)), becomes 

sin 2 9 { ™- shell) = 0.2297 ±0.0019. 

Hence, if we use the experimental measurement of the strange sea asymmetry, the dis- 
crepancy with the standard model is increased to 3.7a significance. 

3.2.3 Nuclear Shadowing 

A recent comment in the literature ^0) has claimed, correctly, that if shadowing were 
significantly different between charged and neutral current neutrino scattering, this would 
affect the NuTeV sin 2 9 W analysis. The authors offer a Vector Meson Dominance (VMD) 
model of shadowing in which such an effect might arise. This model predicts a large 
enhancement of shadowing at low Q 2 which is not observed in deep inelastic scattering 
data. The most precise data that overlaps the low x and Q 2 kinematic region of NuTeV 
comes from NMC which observed a logarthmic Q 2 dependence of the shadowing 
effect as predicted by perturbative QCD. 

Furthermore, shadowing, a low x phenomenon, largely affects the sea quark 
distributions which are common between v and V cross-sections, and therefore cancel 
in R . The NuTeV analysis, which uses v and V data at < Q 2 > of 25 and 16 GeV 2 , 
respectively, is far away from the VMD regime, and therefore the effect even of this 
VMD model is small, increasing the prediction for the NuTeV measured R v and R v by 
0.6% and 1.2%, respectively. Finally, the NuTeV sin 2 9 W data itself disfavors this model 
through its separate measurements of R v and R" ', which are both below predictions, while 
this model increases those very predictions. 

3.3 New Physics 

The primary motivation for embarking on the NuTeV measurement was the possibility 
of observing hints of new physics in a precise measurement of neutrino-nucleon scatter- 
ing. NuTeV is well suited as a probe of non-standard physics for two reasons: first, the 
precision of the measurement is a significant improvement, most noticeably in systematic 
uncertainties, over previous measurements ^> ^> ^\ and second NuTeV's measurement 
has unique sensitivity to new processes when compared to other precision data. In partic- 
ular, NuTeV probes weak processes far off-shell, and thus is sensitive to other tree level 
processes involving exchanges of heavy particles. Also, the initial state particle is a neu- 
trino, and neutrino couplings are the most poorly constrained by the Z° pole data, since 
they are primarily accessed via the measurement of the Z invisible width. 



In considering models of new physics, the "model-independent" coupling mea- 
surement discussed in Section 2 is the best guide for evaluating non-standard contributions 
to the NuTeV measurements. An interesting thing to note about these measurements is 
that they suggest a large deviation in the left-handed chiral coupling to the target quarks, 
while the right-handed coupling is as expected. Such a pattern of changes in couplings is 
consistent with either a hypothesis of loop corrections that effect the weak process itself 
or another tree level contribution that contributes primarily to the left-handed coupling 5 . 
Chiral coupling deviations are often parameterized in terms of the mass scale for a unit- 
coupling "contact interaction" in analogy with the Fermi effectively theory of low-energy 
weak interactions. Assuming a contact interaction described by a Lagrangian of the form 

- C = 2 / - X 2 x YlYIlW.I^h, + IlI^IlW.I^h, +C.C.) , 

H q &{L,R} (A Lfl J 

the NuTeV result can be explained by an interaction with mass scale Aj jL ^ 4 ± 0.8 TeV. 
3.3.1 In teractions from Ex tra U(l) 

Phenomenologically, an extra U(l) gauge group which gives rise to interactions mediated 
by a heavy Z' boson, m Z ' 3> m z , is an attractive model for new physics. In general, 
the couplings associated with this new interaction are arbitrary, although specific models 
in which a new U(l) arises may provide predictions or ranges of predictions for these 
couplings. An example of such a model is an E(6) gauge group, which encompasses 
the SU(3) x SU(2) x £7(1) of the standard model and also predicts several additional 
U(l) subgroups which lead to observable interactions mediated by Z' bosons ^ ' ^) . 
Before the NuTeV measurement, several authors had suggested in the literature that the 
other precision electroweak data favored the possibility of a Z' boson ^5, 36, 37) 

We have analyzed the effect of E (6) -predicted Z' bosons on the NuTeV mea- 
surement of the chiral couplings. As is illustrated in Figure 3, the effect of these bosons 
in the case where the standard model Z and Z' do not mix is primarily on the right- 
handed coupling. It is possible to reduce the left-handed coupling somewhat by allowing 
Z — Z' mixing; however, this possibility is severely constrained by precision data at the 
Z° pole 36 ). 

More generally, a Z' with couplings of the same magnitude as the Z but leading 
to a destructive interference with the Z exchange could explain the NuTeV measurement 

5 It has been noted many times in the literature that the A°p b B deviation, combined with other constraints 
on b quark neutral current couplings implies a shift in the small right-handed coupling. Such a shift is not 
consistent with the hypothesis of a loop-induced correction, either standard or non-standard. 
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Figure 3: The effect of E(6) Z' bosons on the NuTeV measurement of (g^) 2 and (gjf ) 2 . 
The parameter (3 chooses which of the possible U{\) subgroups contributed to the ob- 
served Z'. The standard model prediction is the green point, surrounded by a grid of±la 
top and Higgs mass variations. The upright dark ellipse around shows the effect of an 
unmixed Z'; the lighter ellipse shows the effect of Z — Z' a mixing of 0.003, which is 
already severely constrained by the Z° pole data. 
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Figure 4: Measurements of the neutrino current coupling, interpreted as a neutrino neu- 
tral current interaction rate (oc p^). The precise measurements, T(Z — ► vV) at LEP I 
and the NuTeV measurement of pi are both below expectation. 

if the Z' mass were in the range ~ 1-1.5 TeV. Current limits from the TeVatron experi- 
ments on such Z' are approximately 0.7 TeV 39) Several authors have also recently 
discussed other U(l) extensions in the context of the NuTeV result and found significant 
effects 23 ' 40 ). 

3.3.2 Anomalous Neutrino Neutral Current 

There are few precision measurements of neutrino neutral current interactions. Measure- 
ments of neutrino-electron scattering from the CHARM II experiment ^ and the direct 
measurement of T(Z — > vV) from the observation of Z — > uVj at the Z° pole ^) provide 
measurements of a few percent precision. The two most precise measurements come from 
the inferred Z invisible width ^ ^ ) and the NuTeV result when interpreted as a measure- 
ment of Pq (see Section 2). As is shown in Figure 4, both of the precise rate measurements 
are significantly below the expectation. Although this is not a model-independent obser- 
vation, it is nevertheless interesting to note this connection between two of the discrepant 
pieces of precision electroweak data. 

4 Summary 



The NuTeV experiment has performed a measurement of sin 2 9 W , and finds a deviation 
of three standard deviations from the null hypothesis which assumes the validity of the 



standard model of electroweak interactions. Motivated by the significance of this dis- 
crepancy, we study both conventional and new physics explanations. Several possibilities 
exist, although none is theoretically compelling or has sufficient independent supporting 
evidence to be a clear favorite. Therefore, the cause of this result remains a puzzle. 
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